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SUMMARY 


As  a  sequel  to  the  publication  of  a  preliminary  study  of  warnead 
characteristics  it  was  required  to  modify  certain  of  the  theoretical 
assumptions  made  therein  to  represent  more  precisely  combat  conditions. 
In  this  report,  therefore,  a  theory  is  presented  relating,  generally, 
the  warhead  shape  to  the  distributions  of  fragments  and  targets  about 
the  missile:  it  is  shown  that  the  problem  of  designing  a  warhead  to 
satisfy  a  given  fragment  distribution  is  not  capable  of  a  unique  solu¬ 
tion.  As  an  example  of  the  use  of  the  theory,  and  to  indicate  the 
nature  of  the  results,  it  is  applied  to  a  simple  probability  model 
representative  of  the  engagement  of  a  target  by  missiles  fitted  with 
V.T.  fuses.  Methods  are  suggested  by  which  the  general  theory  may 
be  modified  to  meet  more  stringent  demands. 
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1  Introduction 


An  earlier  report  (Report  of  the  Medium-Range  Anti-Aircraft 
Guided  Y/eapon  Project  Group  -  Appendix’lV)  has  described  an  initial 
and  general  investigation  into  the  desirable  characteristics  of  a 
warhead  for  an  anti-aircraft  guided  weapon  -  the  necessary  warhead 
weight,  charge/case  weight  ratio  and  the  manner  in  which  these  vary 
with  the  accuracy  of  guidance  and  the  conditions  of  engagement.  In 
carrying  out  this  preliminary  study  it  was  convenient  to  make  the 
following  three  assumptions  (as  well  as  others) 

(i)  that  the  performance  of  the  proximity  fuze  was  highly 
idealized; 

(ii)  that  the  dimensions  of  the  target  were  small  compared 
with  othor  distances,  such  as  the  guidance  accuracy, 
involved  in  the  problem;  i.  e#  the  target  could  be 
regarded  as  concentrated  at  a  point; 

and  (iii)  that  the  warhead  could  be  regarded  as  cylindrical  and  yet 

the  angular  distribution  of  its  fragments  could  be  adjusted 
to  match  the  fuze  performance. 

In  the  work  to  be  described  below  it  was  desired  to  modify  these  assump¬ 
tions  so  as  to  be  more  representative  of  actual  conditions  by 

(ij  showing  how  the  actual  burst  pattern  of  any  fuze  could  be 
incorporated  into  calculations,  and  how  the  fragment 
distribution  could  be  matched  to  such  fuze  performance 
so  as  to  give  a  maximum  probability  of  destroying  the 
targe  t ; 

(ii)  introducing  some  parameter  into  the  problem  which  might 
be  used  to  allow  for  the  finite  size  of  the  target; 

and  (iii)  de termining,  approximately,  the  shape  of  the  warhead  case 
which  would  give  the  desired  fragment  distribution  in 
space  -  this  last  being  evaluated  with  the  object  of 
estimating  the  general  form  of  the  warhead,  so  that  its 
installation  into  the  missile  might  be  considered  more 
realistically;  and  so  that,  the  order  of  the  curvature 
of  the  case  being  estimated,  a  study  could  be  initiated 
as  to  the  methods  of  controlling  the  fragment  mass 
distribution  of  such  cases. 

The  burst  pattern  of  the  actual  fuze  which  will  be  used  in  the 
first  designs  of  anti-aircraft  guided  weapons  is  unlikely  to  be  known 
for  some  considerable  time;  this  report  must  therefore  be  regarded  as 
solely  to  establish  a  theoretical  method  of  warhead-fuze  matching,  the 
calculations  presented  here  being  merely  an  example  of  how  such  a  method 
may  be  applied,  an  example  ^/hich  will  be  superseded  when  the  performance 
of  the  finalised  fuze  has  been  determined. 

General  equations  are  established  determining  the  optimum  frag¬ 
ment  distribution,  the  total  number  of  fragments  and  the  probability 
of  destroying  the  target,  all  in  tens  of  the  distribution  of  strikes 
on  vulnerable  targets  about  the  point  of  burst,  which  is  itself 
dependent  on  the  fuze  characteristics  and  the  distribution  of  missiles 
about  the  target  (or  of  targets  about  -the  missile,  the  two  being 
equivalent).  The  equations  are  then  examined  for  a  particular  missile 
distribution  assuming  a  constant  looking-angle  fuze:  the  full  analysis 
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is  shown  in  the  degenerate  case,  corresponding  to  a  looking-angle  of  90°, 
but  in  the  general  case  it  has  been  necessary  to  resort  to  numerical 
methods,  the  resulting  distributions  be  in,,  presented  graphically  to¬ 
gether  wi*i±L  specimen  warhead  contours. 

2  Assumptions 

In  the  general  theory  the  target  is  assumed  to  have  only  one 
component  which  is  singly  vulnerable  and  spherically  symmetrical,  pre¬ 
senting  a  small  area  A:  the  distribution  of  possible  positions  of 
the  vulnerable  area  about  the  point  of  burst  is  supposed  to  be  known. 


To  reduce  the  general  equations  to  a  tractable  form  it  has  been 
necessary  to  construct  a  simple  model  in  the  particular  case.  It  is 
assumed  that  a  fuze  with  constant  looking-angle  y  is  triggered  when 
any  point  D  of  the  target  lies  on  the  surface  of  the  infinite  cone  of 
semi  wangle  Y  about  the  missile  axis  (sec  Fig.  2);  further  the  distance 
from  the  detected  point  to  the  vulnerable  area  (here  supposed  concen¬ 
trated  at  a  point  V)  at  the  instant  of  strike  is  distributed  according 
to  a  throe  dimensional  Gaussian  Error  La w  of  modulus  h  ,  chosen  to 
account  both  for  the  distance  separating  the  detected  and  vulnerable 
points  at  the  instant  of  burst  and,  also,  for  the  motion  of  the  vulner¬ 
able  point  before  the  instant  of  strike.  The  distribution  of  missile 
trajectories  is  a  two  dimensional  Gaussian  one  of  modulus  k  about  a 
central  trajectory  through  0. 


3  Formal  statistical  approach 

It  can  be  shown  that  the  probability  of  fragments  ejected  from 
a  point  striking  an  area  A  at  range  &  is 


1 


-  p  ya2 


(3.01) 


where  p  is  the  fragment  density  in  angular  measure.  Now  suppose  the 
vulnerable  area,  A,  to  be  concentrated  at  a  point  V  whose  co-ordinates 
in  a  spherical  polar  system  referred  to  the  point  of  detonation,  0, 
(Fig.l)  and  the  axis  of  the  missile  uro  (s,  <f>  ,  \(r  ).  Obviously  the 

probability  distribution  of  V,  being  determined  by  the  characteristics 
of  the  fuze  and  the  performance  of  the  guidance  system  will,  in  general, 
not  be  uniform  and,  if  the  probability  of  V  lying  within  the  element 
of  volume  dV  be  expressed  by  dP,  then 


dPq  =  Gi  dV 


(3.02) 


where  is  a  function  of  the  position  in  spaoc  of  V:  expressed  in 

polar  co-ordinates  this  becomes 


aPl  =  gg,  <f>  g  sin  0  ds  <10  a|  (3.03) 

If  the  form  of  g(s,  0  ,  i|r  )  is  determinable  it  is  obviously  desirable 
that  that  warhead  should  be  designed  which  will  relate  the  number  of 
fragments  thrown  in  a  given  direction  to  the  probability  of  V  lying 
in  that  direction:  this  is  equivalent  to  saying 

p  =  *)  (3.04) 

where  the  form  of  f(  0,i|f)  is  to  be  optimised  and  depends  on  g(s,  0  ,  i|r  ). 

A* 
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If  the  area  A  he  singly  vulnerable  the  total  probability  of 
destruction  of  the  target,  P2  is  given  by  the  equation 


2%  % 

W  l  /o(i-a 


s  p 

)  g(s,  i>  ,  t)  sin$  ds  d$  di|r 


r2*  f  f 

a  1  «  /  /  /  e  g(s,^,t)  S  sin0  ds  d <f>  dt 

’  o  '0  1  o 

(3.05) 

The  form  of  f($,  +  )  nay  be  determined  by  maximizing  Po  with  respect  to 


Clearly 

27.  %  ,+  )  *“2 

6P2  =  /  /  /  9  A6f  (0,i|r )  g(s,0,\|r)  sin0  ds  d$  di|r 

o  '  o  ’0 


(3.06) 


to  the  first  order  of  small  quantities. 


Moreover,  the  function  f(0,i|r)  also  determines  the  total  number  of  frag¬ 
ments,  n,  for 


27.  71 

n  =  I  I  f(0,i|r)  sin$  d$  di|r 
o  o 


(3.07) 


Hence,  if  f(0,i|r)  vary  in  such  a  manner  that  n  remains  constant,  it 
follows  that  at  a  stationary  value  of  Pr. 


5Pr  -  a5n  =  0 

where  a  is  an  arbitrary  constant. 
Therefore 


(3.08) 


27.  „7.  M 


•*(*,+)  j? 


0  o 


A8f(^,  \Jr)  g(s,0,  i|/)  sin$  ds  -a8f(^,  \|r)  sin^J  &<f>  d+«  0 


and  since  this  equation  must  hold  for  any  variation  6f  in  f,  it  is 
necessary  that 

~f(M)  4 

As  g(s,0,\|r)<ls  n  a  (3.09) 
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This  is  the  condition  to  ho  satisfied  hy  f(0,i|f)  for  Pq  to  bo  a 
maximum:  it  will  be  seen  that  n  is  dependent  on  a  for  a  defines 


which  in  turn,  defines  n:  therefore,  for  a  fixed  n,  cc  is  also 
fixed.  It  follows  that  if  the  distribution  of  V  about  the  warhead  can 
be  described  by  a  function  E(s,0,\|f)  (which  may  have  to  be  determined 


empirically),  then  the  form  of  f($,ijr)  giving  the  maximum  probability 
of  destruction  may  be  found. 

The  function  so  obtained  defines  the  directions  in  space  in  wnich 
the  fragments  must  travel:  making  allowance  for  the  component  of  velocity 
Imparted  to  the  fragments  by  the  motion  of  the  missile  a  distribution 
function  fs*fcatic  is  derived  representing  the  directions  of  the 

fragment  velocities  relative  to  the  missile. 

It  is  supposed  that  the  direction  of  flight  of  a  fragment  from 
any  point  on  the  missile  surface  is  expressible  in  terms  of  certain 
parameters  depending  on  the  geometry  of  the  warhead,  the  type  of  ex¬ 
plosive  and  the  position  of  the  point  of  initiation.  Conversely  if 
the  fragment  distribution  fs^a^c  ($, \Jr)  is  known  a  v/arhead  contour  may 

be  determined  for  particular  values  of  the  other  parameters:  it  will 
be  shown,  further,  that  such  a  contour  is  not  unique. 

3 . 1  A  particular  case:  the  constant  looking-angle  fuze 

3. 11  The  determination  of 

Suppose  that  the  warhead  is  fused  to  explode  when  any  point  D 
of  the  target,  not  necessarily  a  vulnerable  point,  is  detected  on  the 
surface  of  a  cone  whose  vertex  is  at  0,  axis  coincident  with  that  of 
the  missile  and  semi' -vertical  angle  Y  .  Let  0  be  the  origin  of  a  system 
of  cylindrical  co-ordinates  such  that  D  is  the  point  (2,  r,  0)  and 
V  is  (Z,  R,©)  as  in  Pig,  2. 

•  Then  under  the  assumption  that  the  distribution  of  V  about  0  is 
Gaussian  with  modulus  h, 


and  the  probability  of  V  lying  in  an  element  of  volume  RdRdZd©  when 
D  lies  at  a  given  point  is 


(3.111) 


R  &R  dZ  d.0 


Also  under  the  assumption  that  the  distribution  of  missile  trajec¬ 
tories  is  Gaussian  about  a  central  trajectory  through  D 


Pr.  (  r‘<  PD  <  r  +  dr]  =  2k^e  r  dr 


where  k  is  the  modulus  of  the  distribution.  Hence  the  probability  of 

r  dG  dr 

a  detected  point  lying  \/ithin  the  element  of  area  — — - of  the  surface 


sin  Y 


of  the  cone  is 


r 


(3.112) 


dr  d  3 
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It  follows  that,  if  dP^  is  the  unconditional  chance  of  V  lying  within 
the  element  of  volume  dV,  then,  in  the  notation  of  (3.02) 

-k2r2 

,  ? 

i<re  r 

- ar  ae 

% 


f2%  r  Ir5'-2  ~(h2e  >’^2r2) 

/  --r  G  r  dr  dO  (3.113) 

o  Jo  '^2 


h  = 


-h2e2 


h-’e 


area  oi  cone 


& 


This  function,  expressed  in  polar  co-ordinates  as  g(s,0,\Jf)  will  give 
f(0,i|r)  corresponding  to  a  given  a  by  (3.09). 

N  ow 

l'-  a  (Z  -  z Y  +  (E  cos0  -  r  cos  6 ) "  +  (R  sin  0  -  r  sin  6)^ 
and  2  =  ~ 

V/ 

where  t  —  tan  y 


so  that 


3Pi  - 


2%  «  t  0  -fiit+h2+k2)r2+2h2r  (— +R  cos6cos0+ R  sin0  sin©)  -  h2Z2-h2R2 

1  hrkr  V  * 

e  • 


o  o 


x5/2 


r  d0  dr  1  R  dR  dZ  d0 


x2  co  -(^+h2+k2)r2+ -h2(R2+Z2)  .•^“e  -2h2rR  sinu 

-tf/-*  •  r  •  -J 


f-® 


R  dR  dZ  d© 


% 

where  u  =  —  +0  -0  .  Hence  using  the  relationship 
2 


Mb 


%  +Q  „  -  ^  A 

n  •*  -hj  sin  d 


2  ft 


d6 


-x+c 


may  bo  expressed  in  terms  of  a  Bessel  Function  as 
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t2 


I0(2h2rR)r  dr  j  dR  dZ  d0 


and,  converting  to  polar  eo-ordinatos  vdiere 


=  R  +  Z* 


R  =  s  sin  <f> 
Z  =  S  003  <f> 


dP-|_  = 


f  2lPk2 

l7% 


-(£+h2+I:2)r2+  iR 


COS0  ,2  2 

- h  s 


In(2h2rs  sin$)r  dr 


1 


s 


2 


sin$  ds  d <f>  di|r 


(3.114) 


In  this  case,  accordingly, 


g(s,M)  = 


2h3k2 

~J7T 

%  '  * 


-(>£*: h2+k2)A  lh-r.s.°.°.^-h2s2 


I0(2h  rs  sin$)  r  dr 

(3.115) 


This  function  is,  in  fact,  independent  of  ^ ,  indicating  symmetry  about 
the  missile  axis  of  the  distribution  of  Pi  (this  is  to  bo  expected  owing 
to  the  nature  of  the  fuse).  The  form  of  f(0,i[f)  follows  directly  from 
(3*09)  and  is  found  by  solving  the  equation 


a  ^ 


2h3k2«  r -rmH 


% 


V2 


u 


2 

,v,2  2  2\  2  2h  rs  cos  2  2 

~(j&tr+yr)r  + - r - hV 

4-*—  C 


Ig2h2r  s  sin$)r  dr  J  ds 

(3.116) 


for  given  values  of  a  . 

It  is  easily  shown  that 


VF  sin  ((f)  -<f)~  VM  sin  <f>  (3.117) 

ST3d  C1C 

whore  VM>  Vp  are,  respectively,  the  velocities  of  the  missile  and  of 
the  fragments  relative  to  the  missile; 

<f>  ,  <f>  static  are>  respectively,  the  directions  of  the  velocities  cf 
the  fragments  in  space  and  relative  to  the  missile. 
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Using  this  relation  between  Static  0  it  is  possible  to  derive 

fgtaticC^jt)  from  f(0, i):  the  property  of  symmetry  is,  of  course, 
unaltered. 

% 

3.12  The  solution  when  Y  =  y 


The  integrals  in  equations  3.114,  3.H5  and  3.116  present  some 
difficulty  but  in  the  case  when  Y=y  ,  (that  is,  when  the  fuze  cone 

degenerates  to  a  plane)  the  solution  is  comparatively  simple:  (3.116) 
then  becomes 


a 


2h3k2 


I  -(h^+k2)^ 
/  9 


h2s2 


1 


2. 

I0(2h  rs  sin0)r  dr 


ds 


(3-1201) 


Nov/  it  can  be  shown  that 


I0(at)t  dt  = 


a^ 

e4? 


~(h2+k2)r2 

e  I0  (2n  rs  sin0)r  dr  = - - 


,12 . 2 . 

4h  s  sin  4> 

4(h2+k2) 


2(h2+k2) 


a  = 


h3k2A 


(h^+k2  ) 


.  x  A  .22  h2fs2sin2  <f> 

+  2 - T 


h  +  k 


ds 


Further  it  can  be  shown  that 


Write 


p2  =  h2  fi  - 

v  h2^; 

lY"(h2  cos20  +  k2) 
h2  +  k2 


(3.1202) 
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Then 


<r  =  A  f(<f>, T|r) 


(3.1203) 


h^k2  A 


^/2(h2+k2)  'o 

h\2A 

3/p 

*  2(h2+k2) 


“  _p2s2  _  £ 
a  s2 


ds 


r,2  2 

„  2  P  1 

"  -y - o” 


<3y 


2^v2n  0 


h2k2A 


-2h 


•/a 


\l 


(h2cos2$  +  k2)A  f(^,4 ') 
h2  +.  k2 


v2f 


h^+k2  Jh2oos2$  +k2 


a  = 


2  2 
h  k  A 


-2h. 


2a  J (h2+k2 )  (h2cos2$+k2 ) 


(h^cog^+k2)A  f($,i]r) 

h2  +  k2 


(3.1204) 


It  is  important  to  notice  that  a  limit  is  imposed  on  the  range  of  a  by 
the  modulus  in  the  right  hand  side  of  (3.1202;-).  Therefore 


a  < 


h2k2A 


2  a  y (h2+k2 )  (h2ca  s%  +k2 ) 


(3.1205) 


If  f  h2k2A 

a  <  - 

2A(h2+k2) 

this  condition  is  satisfied  for  all  values  of  but  if 


h2k2A  a  ^  h2kA _ 

2A(h2+k2)  2A./(h2+k2)  (h2.’a  2/  +k2) 


(3.1206) 


it  is  only  satisfied  for  a  certain  range  and  any  integrations  with  respect 
to  '<p  must  he  over  that  range  only,  that  is  between  the  limits  <f>D  and 
%/2  where 

cos  <f>Q  =  i 


hA-jA-A2 


4A2a2 (h2+k2 ) 


-  k 


(3.1207) 
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The  physical  interpretation  of  this  restriction  is,  simply,  that  if  the 
number  of  fragments  is  small  an  optimum  warhead  'would  throw  no  fragments 
in  those  directions  where  a  vulnerable  area  is  rather  unlikely  to  occur, 
but  as  the  number  of  fragments  increases  so  also  does  a  and  the  zone  of 
fragmentation  steadily  widens. 

Rearranging,  (3.1204)  becomes 


,u  ,  Af(0  ^)(h'cos20+k2) 

*J*-h .  I - ; - = 

h2+k2  c 


4ft2a  (h +k2)(h'cos2^+k2) 
iAiAa2 


= 


(h2+k2) 


l6h2A  ( h2c  o  s20  +k2 ) 


loGP 


4^2a2(h2+k2) (h2oos2^+k2)  \  ' 


]AiAa2 


J 


(3.1208) 


or^  in  dimensionless  form, 


f(M)  -  W 

l6a(cos  ^+3")  L 


4 a  (1+3")  (cos"0+3  )\* 


where  the  dimensionless  quantities  3  and  a  are  defined  by  the 
equations 


3  c  |  J  a  =  k2A  . 


Again  the  function  is  independent  of  t  and  is  expressed  here  in  terms 
of  <f>  and  the  parameters  h,  k,  A  and  a. 


It  is  now  possible  to  evaluate  Pp  from  the  equation 


?2 


_A_ 

S2 


<r 


L 


2h\<2 

3/2 


-  (h2 +k2 ) r2-h2  s  2 


I0(2h^rs  sin0)r  dr  J  s^sin^i  ds  d0  d\jr 


(3.1209) 
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Integrating  with  respect  to  r,  as  above, 


a  hks2sin20  2  2 

»  ,  2hV  rz Y  r  a  h2+k2  _hs  2 

P2  “  1  7T  /  J  J  8  - 72T - -  s  sin0  ds  d0 


3/2  . 

ft  3  0  0 


2(h2+k2) 


and  sinoe  f(0,i|r)  is  hero  independent  of  \|r 


P2  =  1  - 


2h\2 


ft  "^2  (h2+k2)  ^0  ^0 


A  h2s2(h2oos2^+k2) 

7  00  -ffo,t)  ^2  ~ - T2T2 - 

J  J  e  +  s2sin0  ds  d$ 


(3.1210) 


Now 


OO 

/ 


2  2  <1 
-ps^ - - 

fit 


s2  s2ds  =  -L 
2P" 


«  P2q2 

-u - JL 

e  u  u2  du 


also 


it  can  be  shown  that,  for  all  values  of  y  , 


— w—  — 

K„  M  =  i  (**)’  /  2_Ji  a„ 


0  W 


v+1 


-p2s2 


_  p  V  1 '  .  _ 

2  s2  ds  =  (jy  K_3y2  (2pq) 


3/2 


P2  =  1  - 


2h3k2 


ft1/2(h2+k2)  VP' 


r  /  *  3/2 
/(* 


K.5^  3in^’ 

(3.1211) 


Now 


K-%  « 


(l + i) 


12. 


d\Jf 
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P2  =  1  -  _ 


2h\2 


V2  1 

^  A  7 C 


%  2  (h“+k2)  0  ^  (4pq)h 


-2pq 


1  +  — —  J  oin$  d^ 


% 


=  1  - 


h3k2 

(h‘  +k2) 


“2pq 

— —  0  sin0  d0  + 


% 


1  ”2P(1  .  J  , 

e  sinp  d <p 


2p3 


] 


Prom  (3.1204)  we  have  the  relation 


h-^k2  “2pq  2%ap 
h2+k2  6  A 


^  % 

2  "  U  ; Sln#  #  + 


% 


2p‘ 


sin  9^  d$ 


(3.1212) 


Now 

% 


—  sin^  d$ 


A  sin^(h^+k^)‘ 

i2r,_„i-j.  , ■ i_2 


(hhk2f 

c  h(h2cos20+k')2  4hA^(h"cos2^+k2) 


2_p?T  a  (A  +k2)  (h^cos^  4-k*~) 


b^kH2 


d^> 


•rc 


(h2+k2 )  /  sin  ^  n  „  ^27ta(h2+k2)  (h2cos2^+k2)'  N 


4h2  (h2cos30+k2) 


2  log. 


h2k2A 


Let 


h  cos  <f>  =  k  tan  0 


Then 

% 


I  -> 


sin0  = 


bm-1 . 1  j. 

2  /  2^a(h2+k2)‘::sec0\ 

4*^  *  tan-1  £  **  °V  h2lcA  /  9 

K 


(h2+k2) 


h2+k2 


?[> 


h3„ 


1  h 
tan  "fe 


logQ  secfi  dfi 


1 
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Also 


% 


x 


2p 


2  sin<f>  d$  = 


(h  +k2)sin$  d$ 
2h2 (h2cos“^>  +k2) 


h2+k2  -1  h 


~  .  27ta(h2+k2) 

sr 2  =  J*  — - -  - - - 

Ah^k 


tan 


-1 


2*a(h2+k2)“ 


1  +  logQ 


h2kA 


tan 


-1 


logg  sec0  d0 


(3.1213) 


This  last  integral  must  "be  evaluated  numerically. 

Again  this  expression  for  P2  must  be  amended  if 


h2k2A 


■<.  a  < 


h2k2A 


2% (h2+k2 )  2%  J (h2+k2 ) (h2cos20+k2 ) 

for  then  (3.C5)  must  be  written 


Pn  a 


2%  x  -<f> 


4 


(l  -  e  s<i)  s"'  sin$  ds  d$  d\|r 


0  J^o  J° 


where  <f>0  is  defined  by  (3.1207)  and,  finally, 


P2  " 


42-2 


hr+k  .  cos  <f> 


Jh2 cos"  ^n+k2 


2fta(ht'+k^)  I,  -1  h  oos$n 
- - - -  ••'tan  - a 

Ah-’k  j  k 


1+  lOgg 


27 


\a  J h2+k2 


h2kA 


r 


-1  hcos$c 


/ 


The  function  /  log  secG  d6  has  been  tabulated. 
o 


log  sec© d0 

(3.1215) 
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3.13  Tho  general  solution 

IThon  Y  is  unrestricted  the  analytical  solution  of  equations 
3.05,  3.114;  3.H5  and  3.116  presents  considerable  difficulty  and  no 
closed  expressions  for  g(s,0,i{r)  f(0,|)  or  P2  have  been  obtained. 

Instead  the  integrals  have  been  evaluated  by  numerical  methods  and 
the  results  are  presented  in  the  form  of  graphs  of  target  and  fragment 
densities  about  the  point  of  burst  together  with  a  -cable  of  P2  against 
a  for  a  particular  value  of  Y  . 

3 . 14  The  barrel  contours 

In  the  special  case  when  i*s-ta-j:ic  ($>1)  is  independent  of  |  the 
warhead  is  a  solid  of  revolutions,  or  barrel  shaped.  Using  the  nota¬ 
tion  of  Pig. 3,  which  illustrates  a  longitudinal  section  through  a  war¬ 
head  initiated  at  I,  it  is  possible  to  equate  two  expressions  for  the 
fragment  density  from  an  elementary  segment  of  thickness  dx,  one  due 
to  the  geometrical  properties  and  one  imposed  by  the  fragment  distri¬ 
bution  function,  namely 


2%y  ds 

4  as 


2%  f 


static 


(0,ijf)  sin0  d <f> 


i.e. 


3y  ds 
4c is  do- 


static 


($,♦)  sin^ 


d£ 

do" 


(3.1U) 


whore  as  is  the  area  of  that  face  of  a  fragment  which  lies  within  the 
warhead  surface  and  on  the  assumption  (previously  made  in  'Report  of 
the  Medium-Range  Anti-Aircraf t  G-uided ’leap on  Project  G-roup'-  Appendix  IV) 
that  75?°  of  the  weight  of  metal  in  the  case  is  converted  to  controlled 
fragments:  s  is  conventionally  measured  countor-cloclcvri.se  round  tho 

contour. 

Prom  the  figure 

4>  =  cr  +0  -  |  (3.142) 


where  <?  is  the  angle  between  the  positive  x-axis  and  the  tangent,  measured 
in  the  positive  direction  of  s,  and  8  is  the  angle  between  the  normal 
to  the  surface  and  the  direction  of  motion  of  tho  fragment  making  tho 
convention  that  8  is  measured  positively  from  the  normal  in  a  counter 
clockwise  direction,  Following  the  work  of  Shapiro  (A  Report  on  Analysis 
of  the  Distribution  of  Perforating  Fragments  for  the  90  mm  M71,  Fused 
T7 4E6,  Bursting  Charge  TNT,  University  of  New  Mexico  UNli/T-234)  the 
following  equation  has  been  adopted  to  define  8 

tan6  =  ~  sec  cr  .  —  (3.143) 

3s  2Vd 

where  V,  Vj)  are  the  initial  fragment  velocity  and  the  velocity  of 
propagation  of  the  explosion  respectively. 

To  maintain  the  contour  equation  in  a  tractable  form  no  attempt 
has  been  made  here  to  allow  for  the  dispersion  of  fragments  about  tho 
direction  ^  although  it  is  hoped  to  examine  the  effect  of  dispersion 
on  the  contour  shape  later. 
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Denoting  — —  by  which  will  be  regarded  as  oonstant,  oquation 


ZVi 


(3.143)  may  be  written 


tan  6  =  j  -2— 


coscr 


(3.144) 


Nov/,  from  the  figure, 


where 


tanP  = 


3  =<r  -  P 


—I -  ,  tancr  =& 

(x  -  x0)  ax 


(3.145) 


(3.146) 


and  the  contour  equation  may  be  derived  by  the  following  method:— 
from  (3.142) 

<10  _  .  d6 

■r  ¥  n  — ■  -L  *T*  * 

dcr  dcr 


where,  by  (3.144)  and  (3.145), 


soc^6  _ £  =  j  seeder  (sin0  +  cos  ersinfl  — ) 
dcr  dcr 


also 


SGC*"0  ~ 
do' 


d0 

cfcf 


2C  d0  dx 
sec  0  —  •  — 

die  dcr 

£5 >2  i  ~ y 

(x-x0)2  +  y2 


2 

seeder 


2 —  using  (3.146) 


dx2 


sin©  sin  3  soc^q- 

_ 

y  -—w- 

dx2 

using  these  results  and  (3.144)  it  follows  that 


by  (3.145) 
and  (3.146) 


66 


(cos2c 7  +  O^COS^p)  L 


3in0  +  sin.0  sin^  seo2cr 


d2y 

dx2 


J 
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from  ( 3 . 141 ) ,  there  fo re 


3y_ 

4&s 


sec'V 

iz 

dx2 


=  f  Static  SW  f1+ - -  (l  4 


4  (cos2c"  +  j^cos^)  V 


y 

dx2 


01 


v/hence,  finally, 


a2y_ 

i?  '  —2. 


coszo'(cos^o*+ j^coo^P+jsinP)  L  4asf s-fca-tic(0  )sin0coso*  Y 


9  9  9 

3y(cos  cr+j  cos  p) 


9 

j  sinQsin  p*'1 


(3.147) 


It  is  important  to  realise  that  the  fragment  distribution  does  not 
define  a  warhead  barrel  uniquely:  in  fact  since  the  differential  equation 
(3.147)  is  of  the  second  order  it  is  formally  possible  to  design  a  doubly 
infinite  set  of  barrels  each  of  which  produces  the  same  distribution  of 
fragments,  the  selection  of  any  one  being  determined  solely  by  installa¬ 
tion  considerations.  The  computed  contour  mil  exist  over  that  range 
of  values  of  <f>  for  v/hich  f  (  0,  i|r)  exists  and  vri.ll  not,  in  general,  termin¬ 
ate  on  the  axis:  if  it  crosses  tho  axis  then  no  real  contour  exists 
satisfying  the  corresponding  boundary  conditions  and  if  it  stops  short 
of  the  axis  the  ends  of  the  barrel  aro  supposed  closed  by  non-fragmenting 
end  plates.  The  boundary  conditions,  which  are  fixed  arbitrarily,  may 
be  considered  as  the  radius  of  the  barrel  and  the  direction  of  the  tangent 
to  its  contour  in  a  longitudinal  section,  both  at  the  point  of  initiation. 

4  Tho  presentation  of  results 

4.1  The  distribution  of  vulnerable  areas  about  the  point  of  burst 


It  has  been  shown  that  the  optimum  fragment  distribution  and  the 
corresponding  Yvarhead  contour  may  be  derived  from  the  function  g(s,0,\|r) 
defining  the  distribution  of  vulnerable  areas  about  the  point  of  burst: 
accordingly  the  first  result  is  a  graphical  representation  of  this 
function  in  the  special  case  of  a  constant  looking-angle  fuse  and  a 
Gaussian  distribution  of  vulnerable  areas  about  the  detected  point 
evaluated  numerically  for  the  following  values  of  the  parameters  (Fig. 4):- 

Y  =  70°  ;  a  =  0.0004  ;  p  =0.25 

which  would  apply  if,  for  example, 

A  =  4  sq.ft. 

R.M.S.  guidance  accuracy  =  100  ft. 

R.M.S.  value  of  l  =  30.6  ft. 

The  distribution  being  symmetrical  about  the  missile  axis,  g(s,0,\|r) 
is  represented  by  contours  in  a  plane  containing  the  axis.  Alternatively 
Pig.  5  displays  the  density  of  vulnerable  areas  as  a  dot  diagram  in  which 
the  number  of  dots,  spaced  conventionally,  in  any  given  part  is  propor¬ 
tional  to  the  probability  of  a  vulnerable  area  lying  v/ithin  the  volume 

17. 
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do  scribed  by  the  rotation  of  that  part  about  the  axis:  the  boundary 
within  which  50?o  of  the  vulnerable  areas  lie  is  also  shown.  It  will 
be  observed  that  this  diagram  shows  reasonable  resemblance  to  the 
burst  positions  actually  obtained  from  V* T.  fuze  trials. 


4 .  2  The  optimum  fragment  distribution  and  probabilities  of  destruc tion 


In  Fig. 6  the  function  f  _.  g  is  plotted  against  0  ,  the 

missile  and  mean  fragment  velocities  being  taken  as  2000  f t/sec  and 
8000ft/sec  respectively.  Corresponding  to  those  values  of  a  values 
of  n  and  Pg  have  been  computed  from  equations  (3.07)  and  (3.05) 
respectively  by  numerical  integration: - 


a 

an  arbitrary  parameter 

the 

V32 

n  * 

total  number  of 
oz.  fragments 

P2 

the  probabilities  of 
destruction  of  the  target 

20  x  10~6 

13300 

0.72 

30  x  10~6 

8700 

0.63 

40  x  10~6 

5900 

0.57 

H  No  significance  should  be  attached  to  the  use  of  l/32  oz.  fragments 
in  this  example  which  is  intended  for  purposes  of  illustration  only: 
at  the  time  the  calculations  wero  performed  it  was  considered  likely 
that  the  optimum  fragment  size  was  approximately  l/32  oz.  ,  but  con¬ 
sideration  is  now  being  given  to  considerably  larger  fragments. 

4. 3  The  barrel  contours 


Having  regard  to  the  complexity  of  equation  (3.147)  it  was  con¬ 
sidered  unreasonable  to  evaluate  more  than  a  few  barrel  contours. 
Accordingly  it  was  decided  to  illustrate  the  variation  in  barrel  shape 
as  the  arbitrary  boundary  conditions  are  changed  by  selecting  a  value 
of  corresponding  to  a  warhead  of  suitable  size  and  then  varying  the 
radius  of  the  barrel  and  the  slope  of  its  contour  at  the  point  of 
initiation:  since  no  other  parameters  wore  varied,  all  the  barrels 

studied  should  produce  identically  the  same  fragment  distribution. 

Figs. 7  and  8  represent  longitudinal  sections  through  two  barrels  which 
are  approximately  at  the  extremes  of  the  possible  range  of  variation. 
Both  weigh  about  150  lbs  and  break  into  8700  fragments  distributed 
according  to  curve  (2)  in  Fig.  6,  but  the  point  of  initiation  lies  at 
opposite  ends  of  the  warhead:  in  each  case  the  constant  j  was  given 
the  value  sin  7°  =  0.122  being  consistent  with  a  deflection  of  7°  from 
the  equatorial  plane,  according  to  the  earlier  theory  of  G.I.  Taylor 
for  the  detonation  of  a  long  cylindrical  cased  charge. 

5  Some  possible  extensions 

It  is  formally  possible  to  modify  the  probability  model  in  a 
number  of  ways  in  order  to  represent  more  closely  combat  conditions. 
Some  of  these  are  briefly  indicated  here  and  it  is  hoped  that  a  more 
detailed  discussion  with  examples  will  bo  published  in  the  near  future. 

No  attempt  has  been  made  in  this  report  to  consider  variations 
in  the  effectiveness  of  fragments  striking  a  vulnerable  area  A:  to  do 
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so  it  is  necessary,  obviously,  to  know/  the  striking  velocity  of  the 
fragment  which  depends  on  a  number  of  variables,  including  the  weapon 
and  target  velocities  and  the  angle  between  them,  the  direction  of 
throw  of  the  fragment,  the  polar  co-ordinates  of  the  target  relative 
to  the  missile,  and  the  height  of  attack.  If  all  these  are  known  the 
proportion  of  striking  fragments  which  do  lethal  damage  may  be  found, 
subject  to  some  assumptions  about  the  distribution  of  their  presented 
areas  at  the  moment  of  impact,  and  a  new  A  defined  as  the  total  vulner¬ 
able  area  reduced  by  that  ratio.  The  quantity  A  is  then  a  function  of 
the  same  variables. 

The  fuze  burst  pattern  may  be  obtained  more  accurately  by  the 
following  method.  On  each  of  a  set  of  polar  diagrams,  representing 
planes  through  the  target  as  pole,  and  for  a  given  direction  of  missile 
attack  in  each  plane,  is  superimposed  a  corresponding  set  of  fuze  burst 
points  (experimental  data)  and  contours  to  show  the  length  of  the  flight 
path  and  necessary  angle  of  throw  for  a  fragment  which  strikes  the  target: 
the  number  of  such  diagrams  must  be  sufficient  to  present  a  fair  picture 
of  the  distribution  of  bursts  in  space  and  associated  with  each  one  there 
wall  bo  the  probability  of  that  direction  of  attack  (dependent  on  the 
method  of  tactical  use  of  the  missile),  and  a  particular  value  of  the 
vulnerable  area  A.  The  ranges  of  flight  and  angles  of  throw,  having 
been  read  from  the  contours,  are  all  transferred  to  one  new  polar  diagram, 
whose  pole  represents  the  position  in  space  of  the  warhead;  the  points 
do  not,  of  course,  carry  equal  weight.  The  new  diagram  is  interpreted 
as  a  density  pattern  such  that  the  number  of  points  falling  within  any 
elementary  area  is  proportional  to  the  probability  that  a  target,  at  the 
instant  a  fragment  reached  the  ring  described  by  the  rotation  of  the 
elementary  area  about  the  missile  axis,  would  lie  within  that  ring.  It 
follows  that  the  diagram  also  represents  the  distribution  function  of  a 
vulnerable  point  about  the  warhead,  denoted  by  g(s,0,\|r)  in  the  particular 
case  of  a  barrel  warhead,  and  that  the  fragment  distribution  function 
f(<£)  may  be  derived  from  it.5* 

YThethcr  or  not  it  is  practicable  to  treat  the  quantity  A  as  a 
function  of  all  the  variables  already  mentioned,  and  yet  to  maintain 
the  problem  in  a  reasonably  tractable  form,  is  still  to  be  determined: 
the  use  of  the  improved  (experimental)  fuze  burst  pattern,  however, 
should  create  no  great  difficulty  and  permits  precise  evaluation  of  the 
manner  in  which  g(s,0,i|r)  depends  on  the  velocities  of  target,  missile 
and  fragment.  To  consider  a  number  of  sub  targets  having  different  vul¬ 
nerable  areas  but  all  situated  at  the  same  point  within  the  target  air¬ 
craft  is  a  simple  extension;  and  it  may  prove  possible  to  allow  for  the 
separation  distances  between  subtargets. 

More  realistic  barrel  contours  may  be  found  by  allovdng  for  the 
dispersion  of  fragments  about  the  predicted  singles  of  throw/:  work  now 
proceeding,  assuming  a  G-aussian  distribution,  will  bo  published  shortly. 
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H  Although  it  is  unlikely  that  consideration  will  be  given  to  the  design 
of  warhead  not  symmetrical  about  the  axis  for  a  number  of  years,  the 
method  may  easily  be  adapted  to  the  general  case,  in  which  f  depends 
both  on  $  and  \|r  . 
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ADDENDUM 

The  comparative  lethalities  of  shaped  and  cylindrical  warheads 


It  was  required  to  estimate  the  advantage  to  be  gained  by  the 
use  of  a  shaped  warhead  instead  of  a  simple  cylinder.  In  this 
addendum  the  lethality  of  a  cylinder  designed  to  contain  the  same 
quantity  of  fragmenting  metal  as  any  one- of  the  family  of  barrel 
warheads,  of  which  those  illustrated  in  Pigs. 7  and  8  are  typical, 
is  calculated  under  the  same  conditions  of  engagement  and  to  the 
some  degree  of  accuracy. 

Method  of  Assessment 

The  area  of  the  curved  surfa.ce  of  the  cylinder  being  fixed  by 
the  requirement  that  both  types  of  warhead  should  produce  the  same 
number  of  fragments  a  length  and  diameter  were  chosen  broadly  similar 
to  those  of  tho  curved  wall  warheads  (but  such  that  the  cylinder  was 
slightly  longer  to  improve  the  spread  of  fragments).  The  character¬ 
istics  of  the  two  typos  are  compared  below,  certain  inconsistencies 
being  the  result  of  a  deliberate  policy  of  favouring  the  cylinder. 

The  warhoad  weights  are  of  the  same  order. 


Barrel1  warhead 

Cylindrical  warhead 

No.  fragments  (l/32  oz) 

8700 

8700 

length 

10  -  12  ins 

17.9  ins 

diameter 

14. 5  -  15.5  ins 

12  ins 

initial  fragment  veloc. 

8000  ft/ soc 

8000  ft/sec 

missile  veloc. 

2000  ft/ sec 

2900  fi/sec 

This  particular  missile  velocity,  in  the  case  of  tho  cylindrical  warhead, 
provides  a  dynamic  fragment  distribution  peaked  at  0  =•  70°  ancl,  according 
to  tho  relations  between  initial  fragment  velocity  and  chargo/case  weight 
ratio  published  in  the  Report  of  the  Anti-Aircraft  Guided  Woapons  Project 
Group,  to  which  reference  has  been  made,  such  a  warhead  would  bo  annular. 


Having  regard  to  the  indefinite  nature  of  annular  warhoad  theory 
at  the  present  time  it  was  decided  to  use  Shapiro's  Formula  to  predict 
the  theoretical  static  fragment  distributions,  denoted  by  f($s),  assuming 
a  solid  charge  initiated  at  the  centre:  in  this  instance  we  have,  in  tho 
notation  of  Pig. 9, 


cot  <f> 


(1) 


whence 


2 

ebe  Rj  cosec  $  ' 

T*  =  ’  (j2  -  cot2^ 


Than,  f($„)  being  definod  as  the  ratio  of  the  number  of  effective  frag¬ 
ments  thrown  from  a  ’sided  of  thickness  dx  to  tho  solid  angle  sub¬ 
tended  at  the  point  of  explosion  by  the  presented  area,  A,  of  the  target, 
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f(0s) 


3j  R  cosec  <f> 

4A  sin  ^(j'bcot2^)'^2 


(2) 


The  geometry  of  the  case  imposes  certain  limits  on  the  range  of  <f>s, 
namely 


jcot  0S|  < 


H 

(e2  +  *¥ 


(3) 


and  it  follows  that  for  a  fixed  value  of  the  product  (t  x  d)  the  frag¬ 
ment  zone  increases  as  Z  increases.  The  function  f(0s)  is  plotted 

in  Pig. 10. 

In  furtherance  of  the  policy  of  favouring  the  cylinder  it  was 
decided  moreover  to  assume  that  fragments  were  subject  to  a  consider^* 
able  Gaussian  error  about  the  predicted  directions  in  a  plane  contain¬ 
ing  the  missile  axis  (the  symmetry  of  the  warhead  made  it  unnecessary 
to  consider  dispersion  in  a  perpendicular  plane);  a  standard  deviation 
of  5°  was  chosen.  The  'dispersed'  static  fragment  distribution  func¬ 
tion  F(0S),  is  compared  with  f($s)  and  with  the  'dispersed'  dynamic 
fragment  distribution  function  P(^(3.)  in  Pigs.  10, 11  respectively. 

It  is  easily  shown  that  the  direction  of  the  fragment  path  in 
space  is  given  by  the  equation 


sin 


0. 364  +  cos  J 


(4) 


and,  hence, 


d-l's  0.3&4-  cos  0d 

d^a  003  (*B  ~  0a) 


(5) 


r 


Since,  obviously,  the  number  of  fragments  in  corresponding  regions  of 
the  static  and  dynamic  distributions  must  be  equal 


P(pE)  sin  <{>s  d$s  =  P(^a)  sin  0d  d0d 
and,  accordingly,  from  equation  (5) 


*(*,)  =  *(*)  aiia  ( x .  °.d5. 

sin  y  cos  (4  -  0a) 

Finally,  under  the  conditions  of  engagement  assumed  in  the 
shaped  warhead  assessment,  namely, 

Y  =  70°;  a  =  0.0004;  P  =  0.25  ; 

the  probability  of  destroying  the  target  represented  by  vulnerable  area 
A  is 
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P2  =  1  - 


4h\2 


oo  oo 


%  0  O^O 


,h2  .2,2  2  2h  rs  cos<fi  2  2 

-(%*h+k)r  + - r - h  s 

t2  J  + 


ID  (2h‘*rs  sin$)r  dr  !  e  s2  s2  sin$  ds  dd 


(7) 


Here  <p  is  written  for  ^  for  brevity. 
Results 


The  fragment  distribution  functions  are  plotted  in  Pigs.  10  and 
11.  Equation  (7)  was  solved  numerically  to  give 


P2  =  0.42 


as  compared  with  a  lethality  of  0.63  for  the  corresponding  shaped  warhead. 
It  is  recognised,  however,  that  the  50fo  increase  in  lethality  obtained 
in  this  particular  case  is  not  necessarily  representative  and  that  in¬ 
spection  of  the  effect  of  varying  other  combat,  target  and  missile  para¬ 
meters,  in  particular  the  fragment  mass,  i3  desirable. 
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FIG.I.  FIG. 2 


FIG. 3. 
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FIG.4. 


constant  looking  anqle  FUZE  (T»70°) 


50*/0  OF  VULNERABLE  AREAS  LIE  WITHIN  THE  HATCHED  ZONE 


THE  FIGURE.  REPRESENTS  THE  GENERATOR  OF  A  THREE  DIMENSIONS  DISTRIBUTION, 
SYMMETRICAL  ABOUT  THE  MISSILE  AXIS  AND,  THEREFORE,  INDEPENDENT  OF  Tjf. 


FIG.4.  CONTOUR  DIAGRAM  TO  SHOW  THE 
DISTRIBUTION  OF  VULNERABLE  AREAS,  g(S.<J>) 
ABOUT  THE  BURST  POINT,  O ,  FOR  A ‘BARREL’ 

WARHEAD. 
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FIG.5. 


<5Wl  p|  2093. 


DOTS  SPACED  CONVENTIONALLY*. SO'lo  OF  VULNERABLE  AREAS  LIE 
WITHIN  THE  BOUNDARY  SHOWN. 


THE  NUMBER  OF  DOTS  IN  ANY  ELEMENTARY  AREA  IS  PROPORTIONAL  TO 
THE  FREQUENCY  OF  VULNERABLE  AREAS  IN  THE  SOLID  OF  REVOLUTION 
GENERATED  BY  THE  ROTATION  OF  THAT  AREA  ABOUT  THE  MISSILE  AXIS. 


FIG.5.  DOT  DIAGRAM  TO  SHOW  THE  DISTRIBUTION 
OF  VULNERABLE  AREAS  ABOUT  THE  BURST  POINT, 

O,  FOR  A ‘BARREL*  WARHEAD. 
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FIG. 6. 


IT*  TO* 
a  *  0*0004 
( 3  » 0*25 

SINCE  THE  WARHEAD  IS  SYMMETRICAL  ABOUT  ITS  AXIS  THE 


FIG. 6.  THE  STATIC  FRAGMENT  DISTRIBUTION  OF 

A 'BARREL'  WARHEAD. 
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FIG.7. 


FIG.7.  BARREL  CONTOUR,  poi  NT  OF  INITIATION  AT  REAR.. 
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FIG. 8.  BARREL  CONTOUR,  point  of  initiation  in  nose. 
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FIG.  9. 


SCALE 


FIG. 9.  THE  CYLINDRICAL  WARHEAD. 
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FIG.IO. 


FIG.IO.  THE  PREDICTED  AND  DISPERSED  STATIC  FRAGMENT 
DISTRIBUTION  FUNCTIONS  FOR  THE  CYLINDRICAL 
WARHEAD  £  @>s)  AND  F  (<$>,) 
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FIG.I  |.  THE  STATIC  AND  DYNAMIC  DISPERSED  FRAGMENT  DISTRIBUTION  FUNCTIONS 

FOR  THE  CYLINDRICAL  WARHEAD,  F(<fc)AND  F(<J>d). 
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